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ABSTRACT 



In the 1975 Hawley and Peebles gave the proposal to use three statistical 
tests for investigations of the galaxies orientation in the large structures. Nowa- 
days, it has been considered as the standard method of searching for galactic 
alignments. In the present paper we analyzed the tests in details and proposed 
a few improvements. Basing on the improvements, the new method of analysis 
of the alignment of galaxies in clusters is proposed. The power of this method 
is demonstrated on the sample of 247 Abell clusters with at least 100 objects in 
each. The distributions of the position angles for galaxies in each cluster are an- 
alyzed using statistical tests: x 2 ? Fourier, autocorrelation and Kolmogorow test. 
The mean value of analyzed statistics is compared with theoretical predictions 
as well as with results obtained from numerical simulations. We performed 1000 
simulations of 247 fictious clusters, each with numbers of galaxies the same as 
the real ones. We found that orientations of galaxies in analyzed clusters are 
not random i.e. that there exists an alignment of galaxies in rich Abell galaxy 
clusters. 



Subject headings: galaxies: clusters: general 
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Introduction 



The analysis of the orien tation o 
galaxies formation scenarios ( 



Doroshkevich 



1973 



Shandarin 



'eebles 



1974 



galaxies' planes is regarded as a standard test o f 



1969 



Dckcl 



Zeldovich 



1985 



1970 



Wesson 



Sunvaew &: Zeldovich 



1982 



Silk 



1983 



1972 



Bower et al. 



20061). Stud ies of the galaxies' planes orientation were conducted as early as in 19th century 



flAbbe 



18751 ) . The review of early methods and re sults of i nvesti gation performed up to the 



Second World War can be found in the article of 



Danverl (ll942f ). today only of historical 



value. The first postwar work, which has been cited to this day, is the treatise of 



Holmberg 



( 119461 ). who compared the numbers of galaxies seen face-on and edge-on, discussed the 
observational effects related to optical measurements of size of galaxy axes, and proved that 
the observed excess of edge-on galaxies is just of observational origin. 

In the early period after the Second World War the researchers were usually 
investigating distributions of position angles within the galax y- rich regions (Cetus, Pisces, 



Hydra, Sex tant, Ursa Maior, Virgo and 



1964 



19681 ). In his two papers Brown (IBrown 



ridanus) of the s 



1964 



sy (jWyatt &: Brown 



1955 



Brown 



19681 ) discovered a departure from 



isotropy in the distributions of posit ion angles . 



angl es of large sem iaxes of galaxies, 
later 



By a naly zing the distributions of pos ition 



Reinhardl (119701 ) and iReinhard fc Roberts! ( 119721 ) and 



Nilsonl (119741 ) found a very weak preference of galaxy plane alignment with the equator 
plane of the Local Supergalaxy. However these results are undermined by the presence of 
background objects. 

Furt her important progress in the investigation of galaxies planes orientation was 



made by 



Hawley &: Peebles 



(119751 ). They discussed in a detailed manner the method of 
investigating the galaxies' orientation through analyzing distribution of position angles as 
well as the influence of possible errors and observational effects. In particular, on the earlier 
papers of Brown, they indicated their insufficient certainty of his results due to possible 
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errors in observations. 



Hawley fc Peebles! (Il975[ ) analyzed the d istributions of posi t ion an gles using x 2 test 



Fourier tests and autocorrelation test. Since 



Hawley fc Peebles! (Il975l) this method w as 



accepted as standard method for ana 



Mac Gilfivrav et al. 



1982 



Diorgovski 



ysis of an galactic alignment (iThomp son 



1983; 



Flin fc Godlowski 



1986 



Kindl 



1987 



1976; 



lin 



1988!: iFfin fc Godlowski! Il989l Il990l Ivan Kampen fc Rheelll990l: ICerne fc Peterson Ill990 



Godlowski! Il993l 119941 iHu et al.lll995l: iGodlowski. Baier fc MacGillivravlll998l: IWu et al 



1997; Godlo wski fc Ostrowski Ill999l: lArval fc Saurer 



Arval. Kandel fc Saurer 



2003; A rval fc Saurer! I2004J; IGodlowski et al.1 12003 lArva 



2006b 



Hu et al 



2006 



200 0l:lBaier. Godlowski fc MacGillivrav 



Wu et al. 



997 



Wu 



fc Saurer 2005a blc. 2006a: 



2006 



Arval. Paudel fc Saurer 



2007, 



Arval 



2008; IGodlowski fc Flinll2010l ; IGodlowski et al.ll2010l : lArval. Bachchan fc Saurer! 12010 



origin al 



1993 



2011) One should n o te th at there ar e several modifications and improvements of 



Hawley fc Peebles! (Il975f) methods 



1994 



Flin fc Godlowsk 



Godlowski. Baier fc MacGillivrav 



1998 



1986 



Arval fc Saurer 



Kindl 



2000 



1987 



Godlowski 



Godlowski et al 



2010 ). The aim o f the p resent paper is to present deeper improvements of the original 



Hawley fc Peebles! (jl975l ) method and show their usefulness for analysis of galactic 
orientations in clusters. The power of this method is shown on the sample of 247 rich Abell 
clusters. 



Following 



Godlowski et al. 



(120051) suggestion tha ; 



richness of the cluster, it was found in 



Godlowski et al. 



alignm ent should increase with 



( I2010l ) that in rich Abell clusters the 



non-randomness of the galaxies' orientation increased with number of objects in clusters. 
The question which arose, is if we could say that in analyzed sample of 247 Abell clusters 
with at least 100 objects each, we found an alignment. For this reason, in the present paper, 
we analyze the distributions of position angles of galaxies belon ging to investigated cluste rs 



using x 2 test, Fourier tests and autocorrelation test applied by 



Hawlev fc Peebles! (!l975h 



(see also 



Flin fc Godlowskil (119861 ): 



Godlowski et al. 



Godlowskil (11993 



1994); 



Godlowski fc FlirJ (12010( 1: 



(120 10l )) as well as Kolmogorow test. For our sample of 247 Abell clusters, 
we compute the mean values of analyzed statistics. Our null hypothesis H is that the 
mean value of the analyzed statistics is as expected in the cases of a random distribution 
of analyzed angles. We compared our results with theoretical predictions as well as with 
results obtained from numerical simulations. 



2. Observational data 



Our observational basis is the same as in iGodlowski et al.l (120101 ) . It is the sample of 



247 A bell clusters with at least 100 objects each, taken from PF catalogue (jPanko &: Flin 



200 61). The structures were extracted from the Muenster Red Sky Survey (MRSS hereafter) 



( jUngruhe et al 



20031 ). MRSS is an optical large scale survey covering the area of 5000 



square degrees in the southern hemisphere with b < —45°. After scanning 217 ESO plates, 
it gives the information about 5,5 million galaxies. The 2D Voronoi tessella tion techniqu e 



was a pplied to the MRSS galaxy catalogue to search for overdense regions (jPanko et al. 



20091 ). PF catalogue, like MRSS, is statistically complete till magnitude value m = 18 m .3 
and it contains structures having at least ten members between magnitude range m 3 and 
m3 + 3 in each structure field. The m 3 is the magnitude of the third brightest galaxy 
located in the considered structure region. The resulting PF includes 6188 such structures. 
We selec t sample of rich c lusters (at least 100 members) being identified with one of ACO 



clusters (1 Abell et al. 



19891 ) . There are 239 such objects in the PF catalogue. Moreover, 



we include 9 objects which can be identified with two ACO clusters, which increase our 



sample to 248 objects. However, we exc 



substructures (IBiviano et al 



1997 . 



ude from our analysis A3822, which potentially has 



20021 ). Therefore, our sample has 247 objects. 



The data for each galaxy member is taken from the MRSS. These includes: the 
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equatorial coordinates of galaxies(a, 5), the diameters of major and minor axes of the 
galaxy image (a and b respectively) and the position angle of the major axis, p. Position 
angles are recomputed from MRSS clockwise system to standard counterclockwise system. 
We perform our comput ation both in Equatorial and Supergalactic Coordinate System 



flFlin fc Godlowski 



1986). Because position angles for the face-on galaxies give only 



marginal information connected with orientation of galaxy, we exclude from analysis all 
galaxies with axial ratio q = b/a > 0.75. 



3. The method of investigation 



The aim of our paper is to check if orientations of galaxies in investigated clusters 
are isotropic. In order to check it we test if the distribution of galaxy position angles 



p (or supergal actic position ang 



introduced by 



H 



Godlowski et al- 



e s P) is isotropic. We apply sta tistical test s originally 



Godlowski! (11993 



19941 ): 



awlev fc Peebles! (119751 ) and later modified by us 

( l2010l ). as well as Kolmogorow test. In all considered tests, the entire 
range of the tested 9 angle (where for 9 one can put p or P respectively) is divided into n 
bins of equal width. In the present paper we use n = 36. 

Let N denote the total number of galaxies in the considered cluster, and N k - the 
number of galaxies with orientations within the fc-th angular bin. Moreover, N Q ^ denotes 
expected number of galaxies in the k-ih bin. In our case all iVo,fe are equal Nq, which is also 
mean number of galaxies per bin. 



Our first test is x 2 test: 



X 



E 



(N k -N Plc ) 
Np k 



E 



(N k - N 0ik ) 



N, 



OA 



(1) 



k=l 1 k=l 

where p k is a probability that chosen galaxy falls into kth bin. We divided entire range 
of a 9 angle into n bins, which gives in the x 2 test (n — 1) degrees of freedom. It means 
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that expected value E(x 2 ) — n — 1 while variance <t 2 (x 2 ) = 2(n — 1). For n = 36 it gives 
E(x 2 ) — 35 and variance cr 2 (x 2 ) = 70. 

We analyzed sample of m = 247 clusters. So we compute the mean value of analyzed 
statistics i.e. in discussed case x 2 value, for whole sample of clusters. If we assume uniform 
distribution of a 9 angle, than expected value E(x 2 ) is again equal 35 while variance 
°" 2 (x 2 ) = ° ~ ' = 0.2834. It gives in our case standard deviation cr(% 2 ) = 0.5324. We check 
our theoretical prediction by numerical simulations. We do this in two ways. 

In the first case we simulate 247 fictious clusters, each with 2360 random oriented 
members galaxies and compute mean value of analyzed statistics (i.e. in the discussed case 
X 2 value). We give 1000 simulations and on this base we obtain: Cumulative Distribution 
Function (CDF) and Probability Density Function (PDF). Expected value of analyzed 
statistics and their variance is computed as well. Comparing them with theoretical 
prediction we are able to check our theoretical assumptions, correctness of the program and 
test quality of used random generator. 

However, please note that the number of galaxies in our real clusters is small in some 



cases, and the x 2 test will not necessarily work well (e.g. t 



re y 2 test requires th e expe cted 



number o: 



' data per bin to equal at least 7; see, however, 



Snedecor fc Cochranl ( 119671 ); 



Domahskil (Il979l ).) For this reason we repeat this procedure but now we simulate 247 
fictious cluster each with number of members galaxies the same as in real clusters. As 
a check, we repeat the derivations for different values of n, but no significant difference 
appear so it is not presented in present paper. 

Now we can compute the mean value of analyzed statistic for real sample of 
analyzed 247 rich Abell clusters and compare it with theoretical predictions and numerical 
simulations. This procedure is also provided for other, presented below, analyzed statistics. 
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The first auto-correlation test quantifies the correlations between galaxy numbers in 
neighboring angle bins. The measure of the correlation is defined as 



C 



Y^{Nk~ N 0tk )(N k+1 - N , k+1 ) 



where N n+1 = N l . 



0,fc+l 



,1/2 



(2) 



Hawley &: Peebles! ( 119751 ) noted that in the case of an isotropic 



distribution, we expect C = with the standard deviation: 

r(C) = n 1 / 2 



cr 



(3) 



Below we show that this result is an approximation which is not valid in our case. 



Hawley fc Peebles! ( 119751 ) result was obtained on the assumption that all N k are 
independent from each other. One should note that E(EX) = S E(X) and moreover, if two 
variables X and Y are independent than we have: E(XY) = E(X)E(Y), D 2 (X + Y) — 
D 2 (X) + D 2 (Y) and D 2 (X Y) = D 2 (X) D 2 (Y) + (E(X)) 2 D 2 (Y) + D 2 (X) (EiY)) 2 . So, if 
all N k are independent than we obtain: 



E(C) ~ E ( V {Nk - N °' k){Nk+1 ~ N °' k+l) ) = E ( ( N ^- N o^)(N k+1 -N 0ik+1 ) 



and 

(W* - iV , fe )(iV fc+1 - JV 0) *+i) \ _ A n2 / (iV fc - iV 0jfe )(iV fe+1 - N 0M 



D 2 (C) = D 2 I ^ ~ ' +1 - °' fc+i; = ^ D 



[N , k N , k+1 ] 1/2 J ^ V [iVA^] 1 / 2 



fc=i 

n /at a t \ /at at \ 1 



= E p2 ( ^ Wt V Wl ) = E 1 = " < 5 > 

fc=l \ iV 0,fc / V iV 0,fc+l / fe=l 

One should note however that the distribution of N k is in fact a polynomial 
distribution and then elements of the covariance matrix of particular N k are given by 
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the formulae Cjj = N Pi(5ij — pj). If two variables X and Y are not independent than 



E(XY) = E(X)E(Y) + Cov(X Y). It leads to a conclusion that now in our formulae for 
E(C) is present an additional term connected with covariance between value of N k and 
N k +i- As a result E(C) = — Yl k =i Np k Pk+i- Because in our case all pk (and as a result 
also Nk,o) are equal {p k = 1/n) than E{C) = — Ylk=i P k = ~Y^k=i = — 1- Moreover 
when we tray to compute D 2 (C) than variance of C contain term which is variance of 
products of Nk and iVfc+i which are not independent. So, correct value of D 2 (C) is different 
from n and is obtained from numerical simulations. 



Differences between our result and lHawley &: Peebles! ( 119751 ) approximation is not 
significant in the case of individual clusters because the difference between results in 
expected value of C (0 or -1) is small with comparison to its standard deviation a(C) ~ \fn. 
However one should note that in our case this difference is important. Our sample has 247 



clusters so standard deviation of C, <r(C) ~ a/^/247 = 0.3818 is significantly smaller than 
a difference in expected values (which is equal to 1). 



If deviat ion from isotropy is a s 



Fourier test ( lHawley fc Peebles 



1975|) 



owly varying function of the angle 9 one can use the 



N k = AV(1 + A n cos20 fc + A 21 sin 29 k ) 



(6) 



We obtain the following expression for the An coefficients 



An 



Y,Z=i(Nk-N , k )cos26 k 
Yl=i N o,k cos 2 29 k 



(7) 



A 



21 



Ek=i(Nk-N 0ik )sm2e k 
J2 n k=1 N 0tk sin 2 29 k 
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These equations are originaly introduced by lHawley &: Peebles! (119751 ). It was written 
in a simple explici t e form in the case when all was equal and and n=36 (equation 25 



Hawlev fc Peeblesl fll975f )l 



Standard deviation of cr(An) and <r(Ai 2 ) is given by expressions: 



ff(A u ) = ^iV , fc cos 2 2^ 



" 1/2 < 2 \l/2 



/ » \ -1/2 2 . J/2 

a(A 21 )= (^iV 0ifc sin 2 2^J = ^— 

The probability that the amplitude 

A 1 =(A? 1 + A 2 1 ) 1/2 
is greater than a certain chosen value is given by the formula 

P(> Ax) = exp (-^JVoAf 
with standard deviation of this amplitude 



1/2 



(9) 



(10) 



ID 



(12) 



(13) 



Th e formula for standard d eviation for a (An), cr(A 12 ), and cr(Ai) (A in the 



original lHawlev fc 



'eeb 



esl ( 119751 ) notation) was also written in a simple explici t form in 



Hawlev fc Peeblesl ( 119751 ) (equation 26). This test was substantially improved by 
( 1994 ) for the case when higher Fourier mode is taken into account: Q 



Godlowski 



N k = N 0tk (l + A n cos 2^ + A 21 sin 20* + A 12 cos40 fe + A 22 sin4^ + ). (14) 

1 However please note that there is a printed error in Godlowski (1994). Eq. 18 should 
have form P(A) = (1 + J/2) exp (- J/2) 
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In our case 



(IGodlowski et al 



all 7 Vo.fr are equal) it leads to formulas for the Ay coefficients 
20101): 



A 



and 



with the standard deviation 



A 2J 



y ~ ELi^ocos 2 2J^' 
££ =1 iV fe sm2J0 fe 



^ =1 iV sin 2 2J^' 



(15) 



(16) 



-1/2 



<t(A]j) = ( ^iV cos 2 2J^ 



,fc=i 



nNn 



1/2 



and 



a(A 2j )=f|:iV sin 2 2^j = (-|_) 



1/2 



(17) 



If we analyze Fourier modes separately, probability that the amplitude 



2 , A2 \V2 



A, = (A? J+ Ay 



is greater than a certain chosen value is given by the formula: 



(19) 



P(> A,) = exp (--iVoA 



(20) 



When we analyze first and second Fourier modes together the probability that the 
amplitude 

A=(A? 1 + A 2 1 + AL + Ay 1/2 (21) 
is greater than a certain chosen value is given by the formulae 



P(>A)=(l + ^ A J 2 )exp(-^iV A 2 ). 



(22) 



The value of coefficient An gives us the direction of departure from isotropy. If, 
An < then the excess of the galaxies with position angles near 90 degrees is observed. 



-12- 



It means that in this case excess of galaxies with the position angel parallel to equatorial 
plane (case p) or parallel to Local Supercluster plane (case P) is observed, while for An > 
the excess of the galaxies with the position angle, respectively perpendicular to equatorial 
plane or Local Supercluster plane is observed. 



In the paper iGodlowski et al.l ( 120101 ) the investigation of the linear regression given 
by y = aN + b counted for various parameters was performed. In the case of position 
angles the linear regression between the values of statistics x 2 ? Ai/cr(Ai), A/cx(A) and 
the number of analyzed galaxies in each particular cluster was studied. It was found that 
non-randomness of galaxy orientation increased with numbers of objects in clusters. To 
test our null hypothesis H that value of the analyzed statistics is as expected in the cases 
of random distribution of analyzed angles, we should now discuss properties of statistics 
Ai/o~(Ai), A/o~(A) as well as properties of the whole Fourier test in more details. 

In one dimensional (ID) case the situation is very clear. Variables An/cr(Aii), 
A2i/cr(A2i) are normalized gausian variables i.e. with expected value equal i.e. 



E(Aij/a(Aij))=0 and Variance equal 1 i.e. D (Ay/cr(A^)) = 1. Of course a (Ay/cr(Ay)) 



1/247 = 0.00405 and a(A ij /a(A ij )) = 0.06363. However, in the case of A 1 / ( x(A 1 ) and 
A/cr(A) variables situation is much more complicated. 



Hawley &: Peebles! (119751 ) (our equation 



In our case, taking into account the equation 26 
T3|) the equations [12] and [20] can be written (in analogy to ID gaussian distribution) in the 
form: 



/ 1 A 2 \ 

P(>A 3 -)=exp(--- 5 ^ y j 



(23) 



while (having in mind that A is given by the equation l2~Tj) the equation [22] could be written 
as: 



P(>A)=(l + i^ |ffl q, 



1 A 2 



2a 2 (A) 



(24) 



Please note however that A,- is described by 2D Gaussian distribution while A is 
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described by 4D Gausian distribution. In an explicit form the equation [12] could be written 
as: 



P(> Ax) = exp 



1 



A 2 



A 2 



(25) 



2 V^ 2 (An) ^(A 21 

A 2 

So, the notation ct2 (a ■> means only that elements of A 2 should be devided by elements of 
covariance matrix Ay. Even more generally it could be written as: 



P(>A 1 )=exp^-i^^G^ 



where I vector is 




(26) 



(27) 



and the matrix G is the inverse matrix to the covariance matrix of Ay (Cov = G 1 ). In the 
4D case vector I has a form: 



A 21 
Ai 2 
y A 22 j 



(28) 



while (using an auxiliary variable J = ^ ^ . JjJj) 

P(> A) = (1 + J/2) exp (- J/2) 



(29) 



Pl ease note t h at th e equation [26] is 2D equivalent of the equation [29] (and the equation 18 



of 



Godlowskil ( 11994] ) paper). Please note that the advantage of our notation is that it could 



be very easily extended on the situation that not all cr(Ajj) are equal and/or non diagonal 
elements of covariance matrix are not disappear (i.e. not all Ay are independent to each 
other). 



One should note that the equation [13] (IHawley &: Peebles! (119751 ) equation 26) is 



obtained as a result of the theorem of propagation of errors. Because in our case 
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<r(Au) = <r(A 2 i) = (2/iV) 1/2 and Ax = (A^ + A 2 2l ) l/2 we obtain the following results: 



dA 



a 2 {A 



lis 



dA 



OA 



21 



<t\A 21 ) 



2A n 



a 2 (A n ) + 



2A 



21 



2^A'i 1 + Al 1 J \2y/A*[+A 
a 2 (A n ) = a 2 (A 21 ) = 2/N 



a 2 (A 



21, 



(30) 



i.e. our equation [131 (or lHawley &: Peebles! (119751 ) equation 26). We obtain analogical 
results in the case when we are taking into account both first and second Fourier modes 
(T 2 (A) = o- 2 (An) = cr 2 (A 21 ) = o- 2 (Ai 2 ) = or 2 (A 22 ) = 2/N. However, please note that the 
theorem of propagation errors is obtained in the linear model. We argue below that a linear 
approximation is not good approximation in our case. 



Please note that a value beeing the sum A?-/cr 2 (Ay) is x 2 distributed (jBrandt 
As a result the value 

A 2 /a 2 (A x ) = A 2 1 /a 2 (A 11 ) + A 2 2l /a 2 (A 21 ) 
has \ 2 distribution with 2 degree of freedom while 

A 2 Ax 2 (A) = A\/a\A xx ) + A 2 1 /a 2 (A 21 ) + A 2 2 /a 2 (A 12 ) + A 2 2 /a 2 (A 22 ) 



19971 ) 



(31) 



(32) 



has x 2 distribution with 4 degree of freedom. So, E(A 2 /a 2 (Ai)) = 2, E(A 2 /a 2 (A)) = 4, 
D 2 (A 2 /cr 2 (Ai) = 4 and D 2 (A 2 /cr 2 (A) = 8. Because we analyzed the A x /(t(Ai) and 
A/cr(A) statistics, the expected value and the standard deviation of this statistics are of 



our interest. The theorem of propagation errors shows that 



a 2 (x) 



dx 



a 2 (x 2 ) 



a 2 (x 2 ) 



a 2 (x 2 ) 



(33) 



dx 1 ) \2\fx 1 J 4x 2 

In our case it leads to results that cr(A 1 /a(A 1 )) = 1/2, and again a (A /a (A)) = 1/2. The 
expected value of E(X) = ^ E(X 2 ) - D 2 (X). In our case 



E 



a(A 1 ) 



E 



Aj 
a 2 (Aj 



D- 



Ai 
cr(A 1 



V2 - 0.5 = 1.2247 



(34) 
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and 



E{ [ A) = r M*>) " fl2 (A. 1 = vCT= 1 - 8708 (35) 



Because in our sample we have 247 clusters than cr 2 (Ai/cr(Ai)) and <7 2 (A/cr(A)) is 

1/2 

equal = 0.002024 while standard deviation of cr(Ai/cr(Ai)) and a (A/ a (A)) is equal 
V 247 = 0-04499. These results are obtained from theorem of propagation errors, so again 
we assume linear approximation. We show in the next section that in presently analyzed 
case it works quite well, but correct values must be obtained from numerical simulations. 

The isotropy of the resultant distributions of the angles 9 can also be investigated using 
Kolmogorov- Smirnov test (K-S test). We assume that the theoretical, random distribution 
contains the same number of objects as the observed one. In such a case statistics A 

A = \/nD n (36) 

is given by limit Kolmogorow distribution, where 

D n = sup\F(x)-S(x)\ (37) 

and F(x) and S(x) are theoretical and observational distributions of 9. Now we can compute 
mean value and standard deviation of analyzed statistic for real sample of analyzed 247 rich 
Abell clusters. Expected value of A, its standard deviation as well as PDF and CDF of A 
are obtained from numerical simulations. 



Numerical Simulations 



A well known problem with random number generators is that their quality is difficult 
to asses in any rigorous way. In fact many of the popular generators used till now failed to 
give correct res ults in multidimensional (sometimes even in two dimensional) simulations 



( ILuescher 



19941 ) . We decided to test a few different number generators to be sure that 
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our result is correct. The first one is built in Fortran Lahey v3 generator. Moreover we 
used classical RANI generator with versio ns (Ranl p t) of Num e rical Recipes ( IPress et al 



19851 ) and more recent version (Ranl nt ) of 



Press fc Teukolskvl fll992h . as well as GGUB S 



subrouti ne from ISML library, minimal standard generator discussed by 



Park & Miller 



(1988((< 



). We compared the result obtained using above generators with a result obtained 



with our basic genera tor which is a subtract-and -borrow random number generator 



originally proposed by 



Marsaglia &: Zamanl ( 



and called RANLUX (level 4) generator (ILuescher 



1991 



) and 



1994 



ater improve d by Martin Luescher 



James 



19941 ) 



The work of Martin Luescher provides the first operational definition of randomness in 
the sense required by Monte Carlo calculations and the Ranlux generator is the first one 
which produces random sequence in which no defect can be observed. The period of the 
generator is about 10 171 . The RANLUX generator is based on a dynamical system which 



may be regarded as a mult i- dimensional version of Arnold's famous cat 
the c at map, the system can be proved to be chaotic in a strong sense 



map. Sim i 



2010). 



Luescher! (11994 , 



arly to 



Fundamentally different between traditional random number generators (TRNG) and 
RANLUX is that TRNGs must be tested because, apart from the testing, there is no reason 
to believe they are at all random. Experience shows that testing is necessary but not 
sufficient. RANLUX, on the other hand, has a good underlying theory, so the purpose of 
testing is only t o make sure that the theory has been understood, applied and programmed 



correctly 



Jamesl (119941 ) . Discussion of diffe rent types of Random G enerator and advantages 



of RANLUX was discussed for example by 



Shchur k Buteral (Il998h . 



We performed 1000 simulations of 247 fictious clusters, each with 2360 random oriented 
members galaxies. As a result we have the sample of I = 1000 values of particular statistics. 
For the present analysis we choose x 2 , C and An/cr(An) statistics for which we have good 
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theoretical predictions only with an exception of the variance for C statistics which is 
obtained from the Peebles approximation. In this Table [TJ we present average value of the 
analyzed statistics, its standard deviation and standard deviation in the sample. Moreover 
we present the stand ard deviation of the standard deviation estimator S which is equal 



a(S) = S/y/2(l - 1) ((Brandt 



19971). 



One should note that RANI generator does not survive the tests. Version of Numerical 



Recipes (IPress et al. 



19851 ) gives wron g result of C statist i cs on more than 40 a level and 
wrong result of x 2 on more than 30 a. IPress &: Teukolskyl (119921 ) version is better but gives 
wrong result of An/o"(An) statistics on more than 10 a level, and more than 2a deviation 
for x 2 statistic. With built in Fortran Lahey v3 generator and GGUBS subroutine, a 
situation is much better, however some deviations from expected value (up to 2a level) are 
observed. The RANLUX generator satisfies all three tests and we choose this generator as 
our base generator. 

In our further analysis we choose six tests. We analyze x 2 > Ai/er(Ai), A/cr(A), C, X 
and An/cr(An) statistics. In the Table |2] we present (as in the Table H]) average value of 
the analyzed statistics, its standard deviation, standard deviation in the sample as well as 
its standard deviation. Because of small number of galaxies in some cluster we repeat our 
analysis (Table |3]) with 1000 simulations of 247 fictious clusters, each cluster with number 
of members galaxies the same as in the real cluster. On the base of these simulations we 
built PDF and CDF presented in the Figures [1] and [2j 

In our procedure we compute the mean values of the analyzed statistics. When the 
errors are normal (Gaussian), what is true at least in the case of A 11 /cr(A 11 ) statistic, that 
parameters are estimated by maximum-likelihood method. They should have asymptotic 
normal (Gausian) distribution. Now we check this suppositions using K-S test. For the test 



we choose the statistics x 2 an d An/cr(An) because we have well theoretical predictions 
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about both expected values and variances of these statistics (one should note however that 
X 2 statistics is \ 2 not normal distributed). As it was shown in the previous sections for x 2 



statistic E(x 2 ) = 35 and a 2 (x 2 ) = 0.2834 while for Au/(r(A u ) statistic E(An/<r(An)) = 



and cr 2 (An/cr(An)) = 0.00405. In order to reject the H hypothesis, that distribution is 
Gaussian with expected value and variance as above, the value of observed statistics A 
should be greater than A cr . At the significance level a = 0.05 the value X cr = 1.358. 

In the case when we performed 1000 simulations of 247 fictious clusters, each with 
2360 random oriented members galaxies (Table [2]) we obtained the values of statistic A 



equal 0.5443 in the case of x 2 statistic and 0.7963 in the case of An /a (An) statistic. 
When we repeated our analysis with 1000 simulations of 247 fictious clusters, each cluster 
with number of members galaxies the same as in the real cluster (Table El) and Figures [T] 
and EJ), we obtained values A = 0.9591 and A = 0.7229 respectively. All these valuese of 
A are significantly less than A cr = 1.358. So we can not exclude our H hypothesis. One 



should note however, that values of x 2 an d A 11 /a(A 11 ) statistics obtained from numerical 
simulations ((Tables El and E]) are a litle bit different from theoretical one. 



Lillieforsl ( 119671 ) showed that the standard tables used for the Kolmogorov-Smirnov 
test are valid when testing whether a set of observations are from a completely specified 
continuous distribution. When we check if the distribution is normal but one or more 
parameters are estimated from the sample then the Kolmogorov-Smirnov test no longer 
applies. At lea st it is not allo wed to use the commonly tabulated critical points. It 



is suggested by 



Masseyl (ll95ll ) that if the test is used in this case, the results will be 



conservative in the sense that the probability of a type I error w ill be smaller than as given 



Lillieforsl ( 119671 ) showed that the results of 



by tables of the Kolmogorov-Smirnov statistic, 
this procedure will indeed be extremely conservative and computed a new table for critical 
value of D — max\F*(x) — Sn(x)\ statistic, where Sn(x) is the observational cumulative 
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distribution function and F*(x) is the cumulative normal distribution function with mean 
value and varianc e estima t ed fro m the sample. His critical values are about 30% less than 
those obtained by iMasseyl (Il95ll ) for classical Kolmogorov-Smirnov test. 



Above modification of the Kolmogorow test is usualy known as Kolmogorow - Lilliefors 
test. For this test the critical value of D cr , at the significance level a = 0.05 for n = 1000 is 
equal 0.028. For analyzed statistics we obtained folowing values of D. For the sample of 247 
fictious clusters, each with 2360 random oriented members galaxies (Table [2]) we obtained 



D = 0.014 in the case of x 2 statistic and D = 0.018 in the case of Au/er(Aii) statistic. For 
sample of 247 fictious clusters, each cluster with number of members galaxies the same as in 
the real cluster (Table |3]) we obtained D = 0.018 in the case of x 2 statistic and D = 0.0135 



in the case of An/o"(An) statistic. Again all these values of A are significantly less than 
critical value and again we can not exclude our Hq hypothesis. 

As a result we can conclude that analyzed statistics can be well described by the 
normal distribution. In particular it means that fluctuations observed in the Figure [2] in 
PDF of x 2 statistics are not in conflict with our prediction that the statistics is normally 
distributed with parameters as obtained from theoretical predictions. 

Independently we analyze situation (Table Hj) when only galaxies brighter than 7713 + 3 
are taken into account. Comparison of the Tables [2], [3] and |4] shows that there are some 
but not big differences between particular cases. When we compare average values of the 
statistics in the Tables [2] and [3] we found that in all cases diferences beetwen them are less 
than theirs 3 standard deviations, however for Ai/ct(Ai), C and A these diferences are on 
the 2a(x) level. It confirmed our sugestion that PDF and CDF used in further analyzis 
should be built on the base of 1000 simulations of 247 fictious clusters, each cluster with 
number of members galaxies the same as in the real cluster. In the case of Tables |3] and H] 
diferences between average values are on 2a{x) level only in the case of A/cr(A) test. So it 
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is not necessary to build special PDF and CDF for this case. 

5. Results 

We analyzed the distribution of the position angles in the sample (A) of 247 rich Abell 
clusters both in Equatorial and Supergalactic coordinate system. Moreover we analyzed 
restricted sample (B) in which only galaxies brighter than + 3 are taken into account. 
The results are presented in the Table EJ 

Our null hypothesis Hq is that mean value of the analyzed statistics is as expected 
in the cases of a random distribution of the position angles, against H% hypothesis that 
analyzed values are different than in the case of random distribution. For the \ 2 test the 
result is significant on 3 a level, for autoccorelation test it is significant on 4 a level, while 
for Fourier test (Ai/er(Ai) and A /a (A) statistics) and Kolmogorow test (A statistics) the 
results are significant on more than 5 a level. In all cases there are no significant differences 
when we analyzed distribution of Equatorial position angles p and Supergalactic position 
angles P. One can see from PDF and CDF presented in the Figures [1] and |2] that probability 
that such results are coming from random distributions is (in all cases) less than 0.1%. 

An/cr(Aii) test does not show any difference from predictions of our null hypothesis 
Hq which means that the average value of the analyzed statistic is as expected in the 
cases of a random distribution. This result, together with the fact that we show no 
significant differences with analysis in Equatorial and Supergalactic coordinate systems, 
shows that observed alignment is not connected with equatorial plane (as expected) nor 
with Supergalactic plane. Interpretation of this conclusion in the context of evolution of 
galaxies in the cluster needs detailed future investigation. 

In our opinion it is because of the influence of environmental effects to the origin of 
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galaxy angular momenta. iGodlowski et all (j201lf ) studied the galaxy alignment in the 
sample of very rich Abell clusters located in and outside superclusters. Even though that 
orientations of galaxies in analyzed clusters are not random, both in the case when we 
analyzed whole sample of the clusters and only clusters belonging to the superclusters, the 
statistically significant difference among investigated samples was found. In contrast to 



whole sample of cluster, where alignment increases with the cluster richness IGodlowski et al. 



( I2010l ) the cluster belonging to the superclusters does not show this effect. Moreover, the 
alignment decreases with the supercluster richness. The observed trend, dependence of 
galaxy alignment on both cluster location and supercluster richness clearly supports the 
influence of environmental effects to the origin of galaxy angular momenta. 

Another im portant possibility is the influence of of the large scale orientation of 



galaxy clusters ( 


HoDkins et al. < 


>005 




Nans et al. 


2009; 


Godlowski & '. 


"lin 


2010: 


Paz et al. 


2011: 


Vera-Ciro et al. 


2011 


Hao et al. 


2011 


. > 


Schaef 


er & Merkel 


2011 




Smargon et al. 


2011 


Blazek. McQuinn & Seliak 


2011 


; Song & Lee 


2011; 


No 


1 & Cohn 


2011 


L) a 


nalized both 
itation of galaxy 


theorticalv and observationalv. IGodlowski & Flin 


( 


20K 


)) studied the oriei 



groups in the Local Supercluster (LSC). It is strongly correlated with t he distribution 



of neighbouring groups in the scale till about 20 Mpc. 



Paz et al. 



(120 111 ) found a strong 



alignment between the projected maj or axis of group shape s and the surrounding galaxy 



distribution up to scales of 30Mpc/h. 



Smargon et al 



(1201 ll ) search for two types of cluster 



alignments using pairs of clusters: the alignment between the projected major axes of the 
clusters founding weak effect up to 20Mpc/h, and the alignment between one cluster major 
axis and the line connecting it to the other cluster in the pair founding strong alignment on 
scales up to lOOMpc/h. 

The change of alignment wi th the surrounding neigbourhood was observed also in 



alignment study in void vicinity (jVarela et al 



201 ll ) being continuation of earlier study of 
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galaxy orientation in regions su rrounding bubble- 



interesting result was found by 



Jones et al. 



ike voids (jTrujillo et al.l 120061 ). Another 



( I2010l ) who reported that the spins of spiral 



galaxies located within cosmic web filaments tend to be aligned along the larger axis of the 
filament, which is interpreted by the authors as "fossil" evidence indicating that the action 
of large scale tidal torques effected the alignments of galaxies located in cosmic filaments. 

For the sample B the results are weaker but still significant. As above, An/cr(An) 
does not show any difference from predictions of our null hypothesis Hq. For x 2 test the 
result is significant at 2a level while for remaining four tests results are significant on more 
than 3<T level. 

Our analysis leads to the conclusion that we observed significant alignment of galaxies 
in our sample of rich Abell clusters. One should note that the most powerful test is the 



t is not a surprise 



jecau se during previous analysis of galactic alignment 



Fourier test 

starting from lHawley &: Peebles! (119751 ) Fourier test was the most sensitive one. Nearly 
the same significance level as Fourier test shows the Kolmogorow test. One should note, 
that in the contrast to analysis of individual structures where the autocorrelation t est 



usually does not lead to a significance conclusion (see for example (IGodlowski 



1993 



1994 



Godlowski et al. 



20051 )). during our analysis the autocorrelation test is more powerful than 



the x 2 test. 



6. Conclusions 



Hawlev k Peebles! (1 19 751 ) for 



We investigated statistical tests originally proposed by 
analysis of the galactic orientations. Basing on analyzed tests, the method of analysis of 
the alignment of galaxies in clusters was proposed. We analyzed the alignment of galaxies 



belonging to 247 Abell clusters containing at least 100 members. The distributions of the 
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position angles for galaxies in each cluster were analyzed using statistical tests: \ 2 test, 
Fourier tests, Autocorrelation test and Kolmogorow test. The mean value of the analyzed 
statistics was compared with theoretical predictions as well as with results obtained from 
numerical simulations. 



The statistical tests originally proposed by lHawley fc Peebles! (119751 ) for analysis of 
the galactic orientations were the % test, the Fourier test and the autoccorelation test. We 
analyzed the autocorrel ation test in more detail a nd some improvements were proposed. It 



was shown that original lHawley fc Peebles! ( 119751 ) result is an approximation which is not 



fully valid in our case. We pointed out that the distribution of the number of galaxies with 
orientations within the fc-th angular bin is in fact a polynomial distribution and then 
in particular Nk are not independent to ea ch other. In the re s ult th e expected value of C 



statistics is equal —1 instead as in original 



Hawley &: Peebles! (1l975[ ) paper. This difference 



is not significant in the case of individual clusters because it is small with comparison to 
its standard deviation a(C) ~ y/n. However in our case when we analyze 247 clusters and 
compute the average value of the statistics this difference begin to be important. This 
is because variance of average values <j(C) ~ a/^/247 = 0.3818 starts to be significantly 
smaller than a difference between our and approximated by Peebles expected values of C . 
Separately we analyzed in detail the Fourier test. We analyzed the properties of whole 
Fourier test as well as the A 1 /cr(A 1 ), A/<r(A) and A 11 /cr(A 11 ) statistics. We compute 
the expected value and the variance of these statistics. The results of our theoretical 
investigations were compared with the numerical simulations. 

Our analysis of the distributions of the position angles of galaxies in rich Abell clusters 
shows that the orientation of galaxies in analyzed cluster is not random i.e. we found an 
existence of alignment of galaxies in the rich Abell galaxy clusters. Five statistical test show 
that distribution of the position angles is not random at least on 3 a level. In all cases there 
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was no significant difference when we analyzed distribution of Equatorial position angles 
p and Supergalactic position angles P. Moreover A 11 /a(A 11 ) statistics do not show any 
significant deviation from randomness. These two facts suggest that observed alignment is 
not connected with equatorial plane (as expected) nor with Supergalactic plane. 



Our previous analysis (IGodlowski fc Flin 



2010 



Godlowski 



2011 



Flin et al. 



20111 ) shows 



the dependency on the alignment of galaxies in clusters and richness of the cluster which 
leads to the conclusion that the angular momentum of the cluster increases with the mass 
of the structure. With such a dependency it is natural to expected that in rich clusters 
significant alignment should be present. In the present paper we confirmed this predictions. 

Usually a dependence between the angu 
is pre sented as empirical relation J ~ M 5 / 3 



ar momentum and t 


re mass of the structure 


Wesson 


1979 


i 


1983 




Carrasco 


1982 




Brosche 



19861 ). In our opinion the observed relation between the richness of the galaxy clus ter 



and the alignment i s due to tidal torque, as suggested by 



Catelan fc Theund (119961 ). Mo reover, the 



Heavens &: Peacock! (119881 ) and 



pointing in the same direction (jNoh &: Lee 



analysi s of the linear tidal torque theory is 



2006a 



bj). They noticed the connection of the 



alignment with the considered scale of the structure. However o ne should note that our 



result is a 



2003, 



so compatible with the prediction of the Li model (ILi 



1998 



Godlowski et al. 



20051 ) in which galaxies form in the rotating universe. 



In our further paper we would like to extend our consideration to analysis of the 
distribution of two angles 5d (the angle between the normal to the galaxy plane and the 
main plane of the coordinate system), and rj (the angle between the projection of this 
normal onto the main plane and the direction towards the zero initial meridian) describing 
the spatial orientation of the galaxy plane. Moreover we would like to investigate if effect 
found in the present paper depends on the cluster BM type. 
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Fig. 1. — The Probability Density Function (PDF) (left panel) and Cumulative Distribution 
Function (CDF) (right panel) for analyzed statistics. The figure was obtained from 1000 
simulations of sample of 247 cluster each with number of members galaxies the same as in 
the real cluster. From up to down we present statistics: x 2 j Ai/cr(Ai), A/a(A). 
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Fig. 2. — The Probability Density Function (PDF) (left panel) and Cumulative Distribution 
Function (CDF) (right panel) for analyzed statistics. The figure was obtained from 1000 
simulations of sample of 247 cluster each with number of members galaxies the same as in 
the real cluster. From up to down we present statistics: C, A, A 11 /cr(A 11 ). 
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Table 1: The comparision of different numerical generators 

Generator Test x a(x) a(x) a{o(x)) 

Lahey X 2 34.974f 0.527f 0.0f66 0.0117 

C -0.9981 0.3743 0.0118 0.0083 

A 11 /(j(A 11 ) -0.0040 0.0616 0.0019 0.0014 

Ranl pt x 2 34.4993 0.5157 0.0163 0.0115 

C -0.5145 0.3569 0.0113 0.0080 

A 11 /(j(A 11 ) 0.0018 0.0613 0.0019 0.0014 

Ranl nr x 2 34.9604 0.5320 0.0168 0.0119 

C -1.0018 0.3836 0.0121 0.0086 

Aii/(t(Ah) 0.0282 0.0632 0.0020 0.0014 

GGUBS x 2 34.9874 0.5505 0.0174 0.0123 

C -1.0110 0.3765 0.0119 0.0084 

An/^An) 0.0021 0.0637 0.0020 0.0014 

RANLUX x 2 34.9978 0.5442 0.0172 0.0121 

C -0.9917 0.3899 0.0123 0.0087 

A 11 /(j(A 11 ) -0.0010 0.0643 0.0020 0.0014 
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Table 2: The result of numerical simulation - sample of 247 cluster each with 2360 galaxies 



Test x a(x) a(x) a(a(x)) 



x 2 


34.9978 


0.5442 


0.0172 


0.0121 


A 1 /a(A 1 ) 


1.2524 


0.0424 


0.0013 


0.0009 


A/a(A) 


1.8794 


0.0460 


0.0014 


0.0010 


C 


-0.9917 


0.3899 


0.0123 


0.0087 


A 


0.7708 


0.0166 


0.0005 


0.0003 


A 11 /cr(A 11 ) 


-0.0010 


0.0643 


0.0020 


0.0014 
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Table 3: The result of numerical simulation - sample of 247 cluster each with number of 
members galaxies the same as in the real cluster. 



Test 


X 


a(x) 


a(x) 


a(a(x)) 


x 2 


34.9798 


0.5364 


0.0170 


0.0120 


Ai/<7(A0 


1.2550 


0.0419 


0.0013 


0.0009 


AMA) 


1.8788 


0.0436 


0.0014 


0.0010 


C 


-1.0195 


0.3749 


0.0119 


0.0084 


A 


0.7720 


0.0168 


0.0005 


0.0004 


An/<7(A n ) 


0.0014 


0.0645 


0.0020 


0.0014 
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Table 4: The result of numerical simulation - sample of 247 cluster each with number of 
members galaxies the same as in the real cluster but only galaxies brighter than m 3 + 3 are 
taken into account. 



Test 


X 


a(x) 


a(x) 


<t(<t(x)) 


x 2 


34.9902 


0.5309 


0.0168 


0.0119 


AiMAx) 


1.2564 


0.0418 


0.0013 


0.0009 


AHA) 


1.8820 


0.0434 


0.0014 


0.0010 


c 


-1.0007 


0.3716 


0.0118 


0.0083 


A 


0.7724 


0.0163 


0.0005 


0.0004 


An/c^An) 


-0.0006 


0.0636 


0.0020 


0.0014 
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Table 5: The value of analyzed statistics real sample of 247 Abell clusters. 



Equatorial coordinates Supergalactic coordinates 



Sample 


Test 


X 


a(x) 


X 


a(x) 


A 


x 2 


36.8591 


0.5924 


36.7899 


0.6315 




A 1 /a(A 1 ) 


1.7046 


0.0622 


1.7021 


0.0626 




A/a(A) 


2.2663 


0.0594 


2.2746 


0.0591 




C 


1.1940 


0.4530 


1.1220 


0.4237 




A 


0.9177 


0.0240 


0.9138 


0.0220 




Aii/<t(Ah) 


-0.0005 


0.0855 


0.0940 


0.0924 


B 


x 2 


36.4000 


0.6072 


36.2919 


0.6124 




A 1 /a(A 1 ) 


1.6283 


0.0577 


1.6316 


0.0578 




A/a(A) 


2.2055 


0.0565 


2.2199 


0.0554 




C 


0.8843 


0.4355 


0.7863 


0.4212 




A 


0.8928 


0.0224 


0.8934 


0.0210 




A 11 /a(A 11 ) 


0.0023 


0.0826 


0.0810 


0.0866 



